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Q_i' Abstract 

(D : 

,^ , We study the production of rj from hadronizing thermal gluons using 

recently proposed rj' — g — g effective vertex. The ry' yield is found to be 
Lj - sensitive to the initial condition. At RHIC and LHC, the enhancement is 

^J ■ large enough to be easily detected. 



1 Introduction 

If the hadronic Lagrangian is symmetric under flavor f/(3) which is spontaneously 
broken, then we would have 9 pseudoscalar Goldstone bosons. In reality, we have 
8 light mesons (vr, K, rf) corresponding to the octet of SU{?)) and one heavy meson, 

V'- 

The pseudoscalar flavor singlet r]' is a remarkable resonance. Its large mass poses 
the Ua{X) problem and its possible resolution relates its mass to the topological 
charge of the QCD vacuum and to the properties of the instanton liquid [|I], 0, |^, Q. 

In heavy ion physics, the r/' meson is a good probe because it has a lifetime 
(1000 fm) long compared to the typical lifetime of a fireball produced by a collision 
of relativistic heavy ions. This was exploited by the authors of Ref. who studied 
possible lowering of the rj' mass by the disappearance of the instanton liquid at high 
temperatures. In it, the authors argued that even in dense matter the t]' meson may 
decouple from the rest of the matter. 

Recently, there was a surge of interest in rj' in the study of S-meson decays and 
search for new physics @, ^, P, |, |10| . In some of these studies, the axial anomaly 



relation 



d^Jl = 2Nj^^Tv[G,.G^^) (1) 

is interpreted to imply that the gluons and rj' have an effective Wess-Zumino-Witten 
type interaction vertex [^ |lOl (See also [p]] ) 

M^,5^' = H{p\ q\ P') 5'^' e,,.p p^ q^ (e^) a (ej), (2) 

where p, q are the gluon momenta and (ep^g)A,7 are the corresponding gluon polariza- 
tion vectors and the superscripts ah denote the color indices of the two gluons. The 
momentum of rj' is denoted by P throughout the paper. By studying J/ip -^ rj' '^ 
decay process, Atwood and Soni |0] found that this process is dominated by on-shell 
gluons and obtained 

Hq = H{0, 0, M^,) ^ 1.8 GeV"^ (3) 

The above ggr]' effective vertex is interesting in many ways. First, since the 
rj' mass is almost 1 GeV, at least one of the gluon momenta involved in the vertex 
should be greater than 0.5 GeV. Therefore the gluon momenta are not soft compared 
to the temperatures achievable in heavy ion collisions. Second, this is a rare occasion 
when we know (at least we can parameterize) how to fuse two on-shell gluons and 



form a hadron. There are models in the hterature that relates constituent quarks to 
the hadrons, but to the author's knowledge, there is no other known matrix element 
between gluons and a known hadron state. 

In this paper, we exploit these unique circumstances and study the production 
of the rj' mesons from a hadronizing quark-gluon plasma. One question we have 
to answer before we proceed is how the interaction strength Hq = if(0,0,M3) 
changes as the temperature increases. In the case of anomalous coupling of photons 
to 7r°, it is known that the coupling strength vanishes in the chiral limit [0, p!3[] 



although the axial anomaly itself is not affected by the temperature |Tj, [T5|, [16 
As the chiral symmetry restoration temperature is approximately the same as the 
deconfinement temperature, one should then ask if the ggrj' strength also vanishes 
as the temperature rises above the critical temperature. 



A partial answer to this question may be given by the result of Ref. |jT3[. In 

^^ 



Ref. [|T^ , the authors carefully analyzed the triangle diagram contribution of n' 
77 decay at finite temperature and obtained 

g-n-y-y = j^e'^g^qqj^lg^gg, a, a In 1/e, ...] (4) 

where mq is the constituent quark mass, e is the electromagnetic coupling constant 
and gTrqq is the 7r° quark anti-quark coupling constant. JF is a finite function of 
coupling constants. If the QCD anomalous coupling of the gluons to rj' is similar to 
7r°77 coupling, the above expression can be rewritten as 

9v'gg = ^0 ~ ^ (5) 

since the coupling constant involved are all 0(1). In a quark-gluon plasma, the u and 
d quark masses vanish. However, the strange quark mass does not vanish. Since 
rj' ~ {uu + dd + ss)/\/3, this indicates that the r]'gg vertex does not necessarily 
vanish in this limit. Rather, it will be proportional to the strange quark mass. 
Furthermore, the r]' mesons produced by the fusing gluons will be dominated by ss 
component. 

There is no doubt that for a more quantitative calculation, we need to calculate 
the rj'gg vertex at finite temperature with full finite temperature complications. In 
this work, we will simply take the coupling to be the same as the vacuum value, 
Ht = 1.8 GeV~^ but take t]' from gluon fusion to be in a ss state so that M^/ fa 
0.7 GeV P] . If we follow the analysis in Ref. , the mass of r]' may be as small as the 
pion mass at high temperature. But in this study, we take the above conservative 
estimate. The finite temperature correction is currently under investigation. 



2 Kinetic Theory Approach 



Kinetic equations are an statement about the change of the phase space density in 
time 

— = (Gain Rate) — (Loss Rate) (6) 

(JjL 

To write down a Boltzmann equation for rj' distribution function, it is easiest to 
start with the decay rate. In terms of the matrix element, the decay rate of rj' to 
two gluons of the opposite colors and different polarizations is given by 

<*"''- - ''i2k I'^^^i^ wk' wk <'''''''^ + ' - ^' <'* 

where p and q are the gluon momenta and P is the rj' momentum. The first factor 
of 1/2 is the symmetry factor. Summing over all final states gives the total decay 
rate. It is then convenient to define 

1^1?'-.,, = E^'^'EI^A.r (8) 

ab A, 7 



It is not hard to show 



using the identities 



|M|;,^^^ = 4|i7orM;^, (9) 

E(s")a(4)a = -^"^ + «pV (10) 

A 

ea/3M.e-'3'"^ = 2(^^<-<^^) (11) 

and the on-shell conditions p^ = q^ = and {p + g)^ = M3. The ap°'p'^ term in 
Eq. (^) does not contribute to M^^ due to the anti-symmetric property of e^p^i,. 

Employing the principle of detailed balance, we then write the Boltzmann equa- 
tion for the phase space density of rj' as 

X |M|^^_^^, [f,{p) f,{q) (1 + f,,{P)) - (1 + f,{p)) (1 + f,{q)) f,,{P)] 

(12) 





(a) 

Figure 1: Feynman graphs for the retarded self-energy of r/'. 

where v = P/Ep. Here it is understood that the distribution functions depend on 
the space time. In the Bohzmann equation, the first term in the colhsion integral 
describes the production of r]' from the gluons and the second term describes the 
decay of rj' into two gluons. These collision terms are essentially the imaginary part 
of the retarded self-energy of rj' [^ depicted in Fig. |I] (b). 



In a series of papers ||T^, Q , it was shown that in a thermal medium, the real 
part of the self-energy must be also included in the mass parameter appearing in 
the Boltzmann equation. With the effective vertex Eq. (0), one can easily calculate 
the one-loop self-energy represented by the Feynman diagrams in Fig.0. The details 
of their evaluation is presented in Appendix ^. In the present case, it turned out 
that the thermal correction is negligibly small up to T ^ 0.5 GeV. Therefore we can 
safely ignore it for our purposes. 

Before proceeding to analyze the Boltzmann equation, we must ask if we can 
use the Boltzmann equation in a quark gluon plasma. In other words, can t]' exist 
in a quark gluon plasma as a quasi-particle? It is possible that an excitation with 
the same quantum numbers as t]' can exist in the plasma (for instance see 0) but 
its width may be too broad to be a quasi particle. To be conservative, we apply the 
Boltzmann equation starting only from one relaxation time before the hadronization 
time unless the hadronization time is shorter than the relaxation time. 

If the steady state is reached during the evolution, then the Boltzmann equation 
dictates that distribution functions become Bose-Einstein functions. In this case, 
the distribution of rj' at the hadronization time will be simply the Bose-Einstein 
distribution with the temperature of T^ ~ 0.17 GeV. This in itself is an interesting 
conclusion because there are surely other hadronization processes and hadronic pro- 



cesses that produce additional tj' . Since r( life time is about 1000 fm and the mean 
free path in dense hadronic matter exceeds 10 fm |0, the final Vj multiplicity will 
exceed thermal model expectation. 

However, local chemical equilibrium between Vj and gluons may not be readily 
reached during the quark-gluon plasma evolution although quarks and gluons do 
reach local equilibrium very fast. Therefore, /^/ must evolve non-trivially in time 
even if the gluons are already locally equilibrated. To calculate such an effect, we 
assume that the gluon density is already thermal and rewrite the above as 

ft/. + V . V/, ^ ji- / p^^ |M|L_, i2.mv ^,-P) 

y-fMUi) li-fAP)/fBE{P)] (13) 

where fsEiP) = l/{e^^^'^ — 1). Substituting the matrix element yields 

dtU + y^-vfr,' = ^^'l^^^' [i-U{P)/fBE{P)]r2{P) (14) 

where the 2-body thermal phase space factor is given by 

The evaluation of r2(P) can be found in Appendix |B[ In the Boltzmann limit, 

T,{P) = ^e-^-/^ (16) 

Svr 

For simplicity, we take the Boltzmann limit from now on. As for the gluon evo- 
lution, we use hydrodynamic models with 1-D expansion to make a simple physical 
estimate. In terms of the space-time rapidity rj = (1/2) ln((t + z)/{t — z)), the flow 
velocity in the 1-D Bjorken model is given by 

u^ = (coshr], 0, 0, sinh?]) (17) 

Using the ideal gas equation of motion e = 3p results in a simple time dependence 
of the temperature 

T{T)=To(^^y (18) 



where r = \/t'^ — z^ is the proper time and Tq is the temperature at the initial 
(proper) time tq. Further, we hmit the Vj momenta to be at the central rapidity so 
that Vz = Pz/Ep = 0. In that case, only the time derivative term from the left hand 
side of the Boltzmann equation remains non-vanishing. The coordinate z and the 
momentum P merely are parameters in 1-D ordinary differential equation 



df{t,z,P) _ \HofM^, 
Jt ~ SttEp 



g-£pcosh„/T(r)_ j^^^^^p^l (^g) 



where we have omitted the subscript label r]' from the distribution and rj and T are 
functions of t and z. This is in the form of a relaxation equation. The relaxation 
time is given by 

where Trei = 4.5 fm is the relaxation time in the rest frame of the rj' and 7p = 
Ep/M^i is the Lorentz 7 factor associated with the rj' momentum. Here we used 
Mri' = 0.7 GeV in accordance with our earlier discussion. Up to the momentum 
of 1 GeV, the typical 7 factor does not exceed 2. Therefore the relaxation time is 
comparable with the typical plasma life time of 1 — 10 fm. The relaxation time tp^ 
is independent of the temperature unless Hq and/or M^i depends strongly on T. 
The solution of the above equation is given by 

f{t,z,P)= |,e-(*-W/o(t',.,P) (21) 

with the initial condition f{Unit,z,P) = and fo{t',z,P) = g-^^^^^^ ,,(*', ^)/T(t',^)_ 
What we are interested in is the distribution function at the hadronization time 
thad- In the Bjorken model, the proper time at the hadronization is given by 

Thad = ro Q-j (22) 

We then take the initial proper time for the rj' evolution to be the larger of tq and 

Tinit ^ Thad ^rel l^oj 

The Mikowskian time and the proper time is related by 



t = Vr2 + ^2 (24) 



Therefore, the farther away from the origin, the later the initial time is. This is due 
to the strong longitudinal flow and time dilation associated with it. The longitudinal 
flow is faster further away from the origin. It also means that at large z, there will 
be very little time between the on-set of r]' production by fusing gluons and the 
hadronization time. 



3 Numerical Results 

To evaluate Eq. (0), we take Tq = 0.334 GeV and tq = 0.6 fm. These parameters are 
taken from a recent hydrodynamic study of the elliptic flow at RHIC[|19[. The initial 
temperature corresponds to the average energy density of about 23 GeV/ fm~^. The 
hadronization proper time with these parameters is Thad = 4.6 fm. Since Thad — ^rei = 
0.1 fm is shorter than tq, we set Tjnit = tq. 

Fig. 1^ shows the numerical solutions with P = within the time interval tjnit < 
t < ^had- In Fig- H we plot the ratio of our solution and what one expects from the 
thermal model, 

/,ac = exp{-Mocoshr]{t,z)/T{t,z)) (25) 

where Mq is the mass of t]' in vacuum. 

The curves in Fig. |^ starts from and keeps growing. This is due to two reasons. 
One, the solution itself overshoots the equilibrium distribution /o (which has the 
same in-medium M^/ = 0.7 GeV) because the temperature is a decreasing function 
of time. Initially the slope df/dt is too steep for the eventual temperature of Tc = 
0.17 GeV. Two, since the in-medium mass is 30% smaller than the vacuum mass, 
/vac in Eq. (|^) decreases much faster than either / or /q as t increases. 

It is also apparent that for larger rj or equivalently larger z, there is not enough 



time between tinit = y ^Init + ^"^ ^nd tgn = JtI^ + ^^ for the solution to grow over 
/vac- Consequently, the enhancement of dN^i /dy at the mid-rapidity {y = 0) 

^ ^' '^ ^2.5 (RHIC) (26) 



{dN,,/dy)^_ 
is not as large as the enhancement of / at r/ = 0. Here we used 



dN^, 



dy 



= f ^ Ep I d'x /(thad, z, Pt) (27) 
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Figure 2: The ratio of the solution of Eq.(|T]) and /o = e-^^ocosh,,{t,2)/T(i,.) at P = 
as a function of time. Here Mq = 0.958 GeV is the vacuum r/' mass. From the 
bottom, the curves correspond to the final space-time rapidities ri = 0.0 through 
0.5 in steps of 0.1 at thad- Calculated with the RHIC parameters. The end points 
corresponds to tinit = jT^iit + ^^ and tgn = Jrl^ + z^ . 



This enhancement factor is not particularly sensitive to the initial temperature. 
Keeping tq fixed, the ratio is 1.8 at Tq = 0.3 GeV, increases up to 2.6 at Tq = 
0.35 GeV and then decreases to 2.2 at Tq = 0.4 GeV. 

One should note that this is on top of other processes that produce rj' at the 
hadronization and in later times. Therefore, this result definitely indicates a large 
enhancement in the 1]' yield due to the thermal gluon fusion process at RHIC. 

At LHC, the initial temperature can reach Tq = 1 GeV. Accordingly, the 
hadronization takes place much later, Xhad = 10 — 20 fm even though the equili- 
bration time is shorter, tq ~ 0.1 fm p2|. Therefore the rate of the change in the 



temperature is slower than the rate at RHIC (recall T = To(ro/r)^/^). This implies 
that even though the initial temperature is much higher than the temperature at 
RHIC, the enhancement factor may not be much different. With tq fixed at 0.1 fm. 



we get 

{dN,./dy)^ 



2-3 (LHC) (2^ 



between Tq = 0.5 GeV and Tq = 1.0 GeV. The maximum enhancement factor 3 is 
reached at Tq = 0.6 GeV. 

At SPS, the enhancement factor is more sensitive to the initial temperature. 



Keeping tq = 0.8 fm ||T9[, the ratio increases as the temperature increases within 



0.2 GeV < To < 0.25 GeV 

{dN„,/dy)f 
0.3 < 'j,r'^ < 1.1 (SPS) (29) 

Since rj' decays to rjTnx in 65 % of the times, one may ask if this SPS result is com- 
patible with the 1] multiplicity measurement by WA80 and the low mass dilepton 
spectrum measured by CERES. Thermal ratio of rj' and r] within the Bjorken sce- 
nario is 17%. Therefore one would expect that about 11 % of 77 comes from rj' decay. 
Doubhng that would indicate about 10% increase in the 77 multiplicity. However, 



at present the experimental uncertainty is bigger than than 10 % IgOj) 21 • 

More detailed information than the yield can be obtained in the transverse mo- 
mentum distribution shown in Fig. ^. There is a clear difference between our calcu- 
lation and the thermal distribution. As one can see, the dependence of the solution 
Eq. (|21| ) on Tq and tq is non-trivial. Since the simple 1-D model we employ does 
not take into account the transverse flow, the slope parameter of the pr spectra in 
Fig. 1^ should be taken as qualitative estimates rather than quantitative predictions. 



4 Discussion and Conclusion 

In this paper, we calculated the yield and the momentum spectrum of the flavor 
singlet rj' mesons produced by the fusion of thermal gluons. It is shown above that 
at RHIC and LHC, there is a significant enhancement in rj' yield. Furthermore, 
the pt spectrum of rj' shows an interesting deviation from the Mt scaling. The on- 
set of the deviation from the naive Mt scaling contains information on the initial 
conditions such as the initial temperature and the thermalization time. This may 
be feasible since rj' has a long life-time and a long mean free path. 

Further implication of our result includes the low mass dilepton enhancement. 
The branching ratio of rj' — > rj-KTi is 65%. A large number of rj' results in a sizable 

10 
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Figure 3: The transverse momentum spectrum of Vj calculated with the solution 
Eq. ( pT] ) (solid line) and the thermal g-^pcoshr^/Tc (^^j^oken line) where E^p is calculate 
with the vacuum r( mass. For the top two lines, we used Tq = 1 GeV and tq = 0.1 fm 
estimated for LHC in Ref.[^. SPS and RHIC parameters are Tq = 0.257 GeV, Tq = 
0.334 GeV and tq = 0.8fm,ro = 0.6 fm respectively ]l|. 
not taken into account. 



Transverse expansion is 



increase in r\ multiplicity which in turn gives rise to the ri Dalitz peak in the dilepton 
invariant mass spectrum. At SPS, the enhancement is not significant enough to 
be noticed. But at RHIC and LHC, there can be a substantial increase of the 
peak. Another observable where 77' enhancement plays a role is the HBT correlation. 
As shown in Ref. ||2^, enhanced r( production reduces the strength of the HBT 



correlation at small p^- 

The above conclusions are based on the effective ggrl vertex deduced by Atwood 
and Soni, the Kinetic (Boltzmann) equation and also the following assumptions: 

(i) The gluons are locally thermalized and follow the hydrodynamic evolution. 

(ii) The strength of ggvl vertex is independent of the temperature. 

(iii) The mass of v( involved in this process is lower than the vacuum value since 
only the s quark loop is involved in the anomalous coupling. 



11 



(iv) Kinetic equation is valid in this regime. 

As to the vahdity of the hydrodynamics, the measurements of the eUiptic flow indi- 
cate the existence of collective motion. Whether this implies thermal and chemical 
equilibrium is not entirely certain. Since our result indicates that the rj' density is 
proportional to the square of the gluon density, the measured t]' dN/dy must reflect 
the underlying gluon distribution be it thermal or the gluon x distribution function. 
For instance, if the plasma is gluon dominated right up to the hadronization, then 
one should see even more t]' than what we have estimated. 

One may also question the validity of the 1-D expansion model we employed. The 
full 3-D calculation with a realistic equation of state is clearly out of the scope of this 
paper. However, our main results should be robust since faster falling temperature 
makes the t]' distribution overshoot even more. 

The temperature dependence of the ggrj' vertex and properties of rj' itself are at 
present not fully understood. In this study, we made simple assumptions to make a 
progress. As we argued in Introduction, the strength of ggt]' vertex will not vanish as 
the strength of tv^'J'J vertex does at Tc. If one accepts the rough estimate Eq. (H) at 
the face value, then the strength of the ggrj' coupling could even be larger than Hq. 
The exact temperature dependence of the vertex, however, has yet to be worked 
out. The crucial question is then: What exactly is the relation between Hq and 
M^/ as a function of temperature? High statistics measurement of rj' spectrum can 
potentially answer this question by employing the ideas developed in this paper. 

The kinetic equation description is valid when the mean free path is much larger 
than any other length scale. This is certainly the case. The mean free time of the 
7]' is longer than 4 fm. The dense medium at T = 300 MeV has much smaller inter- 
particle distance. The effect of inclusion of other processes such as qq -^ r( can be 
roughly estimated by raising the value of ifo- In our calculation, this leads to more 
rl production by shortening trei- 

In summary, we have shown that Vj is a good probe of the gluons density using 
the recently proposed ggr( effective vertex. Other application along the same idea 
includes the investigation of the in-medium properties of Vj and its possible link to 
the fate of the axial anomaly in a quark-gluon plasma. It will be also interesting 
to study formation of v( within gluon jets. These and other aspects are currently 
under investigation. 
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A Real Part of the r]' Self-Energy in Thermal Glu- 
ons 

The Feynman diagrams for the one-loop retarded self-energy of the rj' in equilibrium 
is given in Fig. |I[ These can be calculated in many ways. In this paper, we adopt 



the set of Feynman rules derived in Ref . |2J] . The diagrams then corresponds to the 
expressions 

Sa(P) = Y / (04 (^ (2^)'^(P + ^-P) (Tr M{p, qf) G{p) G{q) (30) 

and 

S.(P) = Y / (01 (01 (2^)'^(P + 1-P) (Tr M(p, qf) A+(p) A+(g) (31) 
where 

G{p) = ^^ + nBE{p)2n5ip') (32) 



14 



and 

A+(p) = ^(p°)27r5(p2) + nB^(pO)27r5(p2) (33) 

The prefactor 1/2 comes from the fact that the two intermediate gluons are identical. 
Using these propagators, it is clear that the real part of the self-energy comes only 
from the diagram (a). Hence we concentrate on evaluation of (a) from now on. 
First, we evaluate the vertex trace 

(TrM(p,g)2) = E'^^'E^a.M,^ 

ab A, 7 

= 16Hi{{p ■ qf - pY) (34) 

using the identities 



and 



EA(epI(4)A = -^"'^ + «pV (35) 

e^p,,e-^P- = 2[gi^g'^^-g^g:;) (36) 

The ap°'p'^ term does not contribute due to the anit-symmetric property of tap^v 
Then 

S,(P) = -S^Hl / ^ ^ (2^)^5(P + q-P)[{p- q? - p\') G{p) G{q) 

= -^^Hl J^,[{P-{P- P)f -P\P- Pf) G{p) G{P - p) (37) 

The zero-temperature part of this diagram is badly divergent. In view of the effective 
theory nature of this vertex, we will simply drop the zero-temperature part in this 
calculation. The thermal part can be separated into two parts, 

S,(P) = S„,i(P) + S„,2(P) (38) 

where 

E„,(P) = -16.i/oV|^(p.i-)' ,,,_,;,^^,^ n(ri2..(p') (39) 

and 

S„,2(P) = -SiHlj^^{p-Pfn{Ep_,)2n5{M^-2P-p)n{p)27^5{p^){AQ) 
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and we used the on-shell conditions p^ = and P^ = M^. 
The real part of the self-energy comes only from Sa,i(P). 

ReS,et(P) = ReS„,i(P) 



«^«/(5^(^'^''pp 



M2 - 2P • 



p 



n{p)2iT6{p^ 



-AH^M^ 



+ AHlM^ 



d'^p 
f d^p 



(27r) 



n{p)2n5{p'^ 



■PP 



M2 - 2P ■ p 



n{p)2Tx5{p'^ 



(41) 

where PP signifies the principal part. 

For simplicity, we orient P^ = {Ep, 0, 0, P) and approximate the Bose-Einstein 
factor by a Boltzmann factor 



n{p) ~ e 



-p/T 



(42) 



Then 

ReS(P) 



-mlM'j^^jdppe-\^\'^ 



+ AHlM* 



d*p 

(2^ 



PP 



M2 - 2P ■ 



=-bl/r 



P 



27r5(p^ 



-4i^o^^^7^,. 



(2vr)^ 



+ 4^/^--^^^{^^ 



M^ - 2£;p + 2Pj9 
M2 - 2Ep - 2Pp 



+ ln 



M^ + 2Ep + 2Pj9 
M2 + 2Ep - 2Pp 



The result can be expressed in terms of the exponential integral funct 

ReSret(P) ^ 
where 



ions 



(27r)2 



M2 

/-4PT^(^) 



(43) 



(44) 



- e^-/^ Ei f ^^— ^l + e^/^ Ei f Z^±^ 
V 2T / V 2T 



+ e(P+i?p)/TEi 



(P + ^p) \ _ J.J /P+^ 
2T / V 2T 



(45) 
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Numerically, between T = 0.17 GeV and T = 1.0 GeV, the second term is important 
only near P = 0. Therefore we approximate the above with 

ReSret(P) ^ Sret(P = 0) 

= -M^^^ (l-^ (e^/2TEi(_^/2T) + e-^/'^Ei{M/2T))) 

(46) 

We can self-consistently determine M by solving the following equation for M with 
Mo = 0.7 GeV: 

M2 = M2 + ReSret(0) (47) 

Numerically solution of this equation indicates that M{T) is a slow varying function 
of T. At T = 0.17 GeV, M{T) = 0.7 GeV is indistinguishable from Mq. Even at 
T = 0.5 GeV, M(T) = 0.68 GeV. Only around T = 1 GeV, M{T) = 0.64 GeV 
is appreciably different from Mq. However, at this temperature, our Boltzmann 
approximation is no longer appropriate. We can ignore the temperature dependence 
of the f]' mass for our estimates. 

B 2 Body Thermal Phase Space 

We start from the expression 

^='^' = / (2^ (202^ P')'^'^ + " - ^'^»<''' ^»''' '«> 

Carrying out the d'^p integral and the q angle integral yields 

^^in = \ I (2^)3^^jp _ q| (2^)'^(|P - ql + g - Ep)mP - ql) fM 

11 1 /•'Jmax 



where 



and 



11 1 /-(Jmax 

I^r^ / dq f.iEp - q) f,{q) (49) 

4 27r F J<i,„,i„ 



M2 
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M2 



Qn 



2{Ep - P) 

Using the Bose-Einstein functions for /, we get 

T 



(51) 



r2(P) = fsEiEp) 



in P 



In ( e«--/^ - 1 ) - In ( e«--/^ - 1 



+ In fe^^/^ - e'^--/^) - In fe^^/^ - e^-^'^/^ 



In T -^ limit, we recover the Boltzmann result 



(52) 



UP) 



-Ep/T 



ilT 



(53) 



